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^ , Microscopic studies of nuclear matter under diverse conditions of density and asymmetry 

Ch ■ are of great contemporary interest. Concerning terrestrial applications, they relate to 

future experimental facilities that will make it possible to study systems with extreme 
neutron-to-proton ratio. In this talk, I will review recent efforts of my group aimed 
at exploring nuclear interactions in the medium through the nuclear equation of state. 
The approach we take is microscopic and relativistic, with the predicted EoS properties 
derived from realistic nucleon-nucloon potentials. I will also discuss work in progress. 
Most recently, we completed a DBHF calculation of the A hyperon binding energy in 
nuclear matter. 
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' 1. Introduction 

. ^ , The properties of hadronic interactions in a dense environment is a problem of 

' fundamental relevance in nuclear physics. Such properties are typically expressed 

^ . in terms of the nuclear equation of state (EoS), a relation between energy and 

density (and possibly other thermodynamic quantities) in infinite nuclear matter. 
The infinite geometry of nuclear matter eliminates surface effects and makes this 
system a clean "laboratory" for developing and testing models of the nuclear force 
in the medium. 

The project I review in this talk is broad-scoped. We have examined several EoS- 
sensitive systems/phenomena on a consistent footing with the purpose of gaining 
a broad overview of various aspects of the EoS. We hope this will help identify 
patterns and common problems. 

Our approach is microscopic, with the starting point being realistic free-space 
interactions. In particular, we apply the Bonn potential in the Dirac-Brueckner- 
Hartree-Fock (DBHF) approach to asymmetric nuclear matter as done earlier by 
the Oslo group. ^ The details of the calculations have been described previously.^ As 
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it has been known for a long time, the DBHF approximation allows a more realistic 
description of the saturation properties of symmetric nuclear matter as compared 
with the conventional Brueckner scheme. The leading relativistic effect characteristic 
of the DBHF model turns out to be a very efficient saturation mechanism. We 
recall that such correction has been shown to simulate a many-body effect (the "Z- 
diagrams") through mixing of positive- and negative-energy Dirac spinors by the 
scalar interaction.'* 

In what follows, I will review some of our recent results and discuss on-going 
work. 1 stress again that these efforts belong to the broader context of learning more 
about the behavior of the nuclear force in the medium using the EoS of infinite 
matter (under diverse conditions of isospin and spin asymmetry) as an exploratory 
tool. I also emphasize the importance of fully exploiting empirical information, 
which is becoming more available through collisions of neutron-rich nuclei. 



2. Isospin-asymmetric nuclear matter 

2.1. Seeking laboratory constraints to the symmetry potential 

Reactions induced by neutron-rich nuclei can probe the isospin dependence of the 
EoS. Through careful analyses of heavy- ion collision (HIC) dynamics one can iden- 
tify observables that are sensitive to the asymmetric part of the EOS. Among those, 
for instance, is the neutron-to-proton ratio in semiperipheral collisions of asymmet- 
ric nuclei at Fermi energies.^ 

In transport models of heavy-ion collisions, particles drift in the presence of 
an average potential while undergoing two-body collisions. Isospin-dependent dy- 
namics is typically included through the symmetry potential and isospin-dependent 
effective cross sections. Effective cross sections (ECS) will not be discussed in this 
talk. However, it is worth mentioning that they play an outstanding role for the 
determination of quantities such as the nucleon mean-free path in nuclear matter, 
the nuclear transparency function, and, eventually, the size of exotic nuclei. 

The symmetry potential is closely related to the single- neutron/proton poten- 
tials in asymmetric matter, which we obtain self-consistently with the effective in- 
teraction. Those are shown in Fig. [T] as a function of the asymmetry parameter 
a = £2Ji£z_ The approximate linear behavior, consistent with a quadratic depen- 
dence on a of the average potential energy, is apparent. Clearly, the isospin splitting 
of the single-particle potential will be effective in separating the collision dynamics 
of neutrons and protons. The symmetry potential is shown in Fig. [2] where it is 
compared with empirical information on the isovector part of the optical potential 
(the early Lane potential^). The decreasing strength with increasing energy is in 
agreement with optical potential data. 
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Fig. 1. The single-neutron (upper panel) and single-proton (lower panel) potentials as a function 
of the asymmetry parameter for fixed average density {kp = 1.4/m~^) and momentum (k = kp). 




Ekin (MeV) 

Fig. 2. The symmetry potential as a function of the nucleon kinetic energy at kp = lAfm~^. 
The shaded area represents empirical information from optical potential data. 

2.2. Summary of EOS results and comparison with recent 
constraints 

In Table [T] we summarize the main properties of the symmetric and asymmetric 
parts of our EoS. Those include saturation energy and density, incompressibility if, 
the skewness parameter K' , the symmetry energy, and the symmetry pressure, L. 
The EoS for symmetric and neutron matter and the symmetry energy are displayed 
in Figs. HH 

A recent analysis to constrain the EoS using compact star phenomenology and 
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Fig. 3. Energy /particle in symmetric (solid) and neutron (dashed) matter, as a function of density 
(in fm~^). 




Fig. 4. Symmetry energy as a function of density (in fm •^). 



HIC data can be found in Ref.^ While the saturation energy is not dramatically 
different between models, the incompressibility values spread over a wider range. 
Major model dependence is found for the K' parameter, where a negative value 
indicates a very stiff EoS at high density. That is the case for models with parameters 
fitted to the properties of finite nuclei, whereas flow data require a soft EoS at the 
higher densities and thus a larger K' . The L parameter also spreads considerably, 
unlike the symmetry energy which tends to be similar in most models. For L, a 
combination of experimental information on neutron skin thickness in nuclei and 
isospin diffusion data sets the constraint 62 MeV < L < 107 MeV. 

Overall, our EoS parameters compare reasonably well with most of those con- 
straints. They also compare well with those from other DBHF calculations reported 
in Ref.^ 

The parameters in Table [T] are defined by an expansion of the energy written as 



with e = [p ~ Pa) / Ps, and ps is the saturation density. 
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Tabic 1. An overview of our 
predicted properties for the 
EOS of symmetric matter and 
neutron matter. 



Parameter 



Predicted Value 




Ps 

K 
K' 



-16.14 MeV 
0.185 fm-3 
259 MeV 
506 MeV 
33.7 MeV 
70.1 MeV 



3. Spin-polarized neutron matter 

In this section we move on to a another issue presently disciisscd in the literature 
with regard to exotic states of nuclear /neutron matter, namely the aspect of spin 
asymmetry and (possible) spin instabilities. 

The problem of spin-polarized neutron/nuclear matter has a long history. Ex- 
tensive work on the this topic has been done by Vidaha, Polls, Ramos, Bombaci, 
Miither, and more (see bibliography of Ref.^ for a more complete list). The ma- 
jor driving force behind these efforts is the search for ferromagnetic instabilities, 
namely the existence of a polarized state with lower energy than the unpolarized, 
which naturally would lead to a spontaneous transition. Presently, conclusions differ 
widely concerning the existence of such transition and its onset density. 

A coupled self-consistency scheme similar to the one described in Rcf.^ was de- 
veloped to calculate the EoS of polarized neutron matter. The details are described 
Ref.^ As done previously for the case of isospin asymmetry, the single-particle po- 
tential (for upward and downward polarized neutrons), is obtained self-consistently 
with the effective interaction. Schematically 



with i,j=u, d. The nearly linear dependence of the single-particle potentials on the 
spin-asymmetry parameter* 



with /? the spin-asymmetry parameter, is reminescent of the analogous case for 
isospin asymmetry and may be suggestive of a possible way to seek constraints 
on Us, the "spin Lane potential", similarly to what was discussed above for the 
isovector optical potential. Namely, one can write, for a nucleus, 




(2) 



(3) 



t/»[/o + C/a(s-I])M, 



(4) 



with s and S the spins of the projectile nucleon and the target nucleus, respectively, 
and extract an obvious relation with the previous equation. (In practice, the situa- 
tion for an actual scattering experiment on a polarized nucleus would require a more 
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Fig. 5. Spin symmetry energy as a function of the neutron density (in fm ^). 

complicated parametrization than the one above, as normally a spin-unsaturated 
nucleus is also isospin-unsaturated.) 

As already implied by the linear dependence of the single-particle potential dis- 
played in Eq. ([3]) , the dependence of the average energy on the asymmetry parameter 
is approximately quadratic, * 

e{p,P)^e{p,Q) + S{p)p\ (5) 

where S{p) is the spin symmetry energy, shown in Fig. [5] The spin symmetry energy 
can be related to the magnetic susceptibility through 

The rise of S{p) with density shows a tendency to slow down, a mechanism that we 
attribute to increased repulsion in large even-singlet waves (which contribute only 
to the energy of the unpolarized state) . This could be interpreted as a precursor of 
spin instability. In Table [2] we show predicted values of the ratio Xf/Xj where xf 
is the susceptibility of a free Fermi gas. 

Concerning the possibility of laboratory constraints which may help shed light 
on these issues, magnetic properties are of course closely related to the strength of 
the effective interaction in the spin-spin channel, which suggests to look into the 
Go Landau parameter. With simple arguments, the latter can be related to the 
susceptibility ratio and the effective mass as 

X ^ 'm*/m 

XF 1 + Go ^ ' 

Thus, Go< -1 signifies spin instability. (Notice the analogy with the formally similar 
relation between the incompressibility ratio K/Kp and the parameter _Fbj where 
Fq< -1 signifies that nuclear matter is unstable.) At this time, no reliable constraints 
on Go are available, due to the fact that spin collective modes have not been observed 
with sufficient strength. 
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Tabic 2. Ratio 
Xf/x s-t differ- 
ent densities. 



p/po 



Xf/x 



0.5383 
0.8886 
1.3651 
1.9873 
2.7741 
3.7457 



2.1455 
2.3022 
2.4950 
2.6411 
2.6824 
2.6548 



In closing this section, we note that we found similar considerations concerning 
the trend of the magnetic susceptibility to hold for symmetric nuclear matter as 
well as for neutron matter. 

4. Work in progress: non-nucleonic degrees of freedom 

There are important motivations for considering strange baryons in nuclear mat- 
ter. The presence of hyperons in stellar matter tends to soften the EoS, with the 
consequence that the predicted neutron star maximum masses become considerably 
smaller. With recent constraints allowing maximum masses larger than previously 
accepted limits, accurate calculations which include strangeness become important 
and timely. On the other hand, remaining within terrestrial nuclear physics, studies 
of hyperon energies in nuclear matter naturally complements our knowledge of finite 
hypernuclei. 

The nucleon and the A potentials in nuclear matter are the solution of a coupled 
self-consistency problem, which reads, schematically 



To confront the simplest possible scenario, one may consider the case of symmetric 
nuclear matter at some Fermi momentum kp in the presence of a "A impurity", 
namely kp w 0. Under these conditions, the problem stated above simplifies con- 
siderably. Such calculation was done in Ref.^ within the Brueckner scheme. 

We have done a similar calculation but made use of the latest nucleon-hyperon 
(NY) potential of Ref.,^" which was provided by the Jiilich group. In a first ap- 
proach, we have taken the single-nucleon potential from a separate calculation of 
symmetric matter. (Notice that the A potential is quite insensitive to the choice 
of Un, as reported in Ref.^ and as we have observed as well.) The parameters of 
the A potential, on the other hand, are calculated self-consistently with the Gnk 
interaction, which is the solution of the Bethe-Goldstone equation with one-boson 
exchange nucleon-hyperon potentials. In the Brueckner calculation, angle-averaged 
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Pauli blocking and dispersive effects are included. Once tlie single-particle potential 
is obtained, the value of —U\{p) at p=0 provides the A binding energy in nuclear 
matter, B\. 

As shown and discussed extensively in Rcf.,"'^'' there arc several remarkable dif- 
ferences between this model and the older NY Jiilich potential, and those seem to 
have a large impact on nuclear matter results. The main new feature of this model 
is a microscopic model of correlated tttt and KK exchange to constrain both the tr 
and p contributions.^" With the new model, we obtain considerably more attrac- 
tion than Reuber et al.,^ about 49 MeV at fc^=1.35 fm^^ for Ba- We have also 
incorporated the DBHF effect in this calculation (which amounts to involving the 
A single-particle Dirac wave function in the self-consistent calculation through the 
effective mass) and find a moderate reduction of i?A by 3-4 MeV. A detailed report 
of this project is forthcoming. 

The natural extension of this preliminary calculation will be a DBHF self- 
consistent calculation of Un, Ua, and U-e for diverse A and S concentrations. 

5. Summary and conclusions 

I have presented a summary of recent results from my group as well as on-going work. 
Our scopes are broad and involve several aspects of nuclear matter, the common 
denominator being the behavior of the nuclear force, including its isospin and spin 
dependence, in the medium. I have stressed the importance of seeking and exploiting 
laboratory constraints. In the future, coherent effort from theory, experiment, and 
observations will be the key to improving our knowledge of nuclear matter and its 
exotic states. 
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